Critical properties of ID spin-gapped fermions at the onset of magnetization 
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We develop an effective field theory for finding critical properties of ID spin gapped fermions at 
the onset of magnetization. It is shown how the spin-charge interaction leads to a linear critical 
behavior and finite susceptibility for a wide range of models. We also discuss possible manifestations 
of spin-charge coupling in cold atomic gases. 



Spin-charge separation is a distinctive feature of one di- 
mensional (ID) multicomponcnt fermionic/bosonic sys- 
tems. In contrast to the Fermi liquid picture, elementary 
excitations in ID systems are not quasiparticles carry- 
ing both spin and charge, but rather spin and charge 
waves that propagate with different velocities [1|. This 
has been addressed in a number of experimental stud- 
ies and demonstrated, in particular, in experiments with 
quantum wires in semiconductors |2j] . Currently, there is 
a growing interest in revealing effects of spin-charge sep- 
aration in experiments with cold Fermi gasesl3[, where 
the ID regime has been recently achieved |j,l5(. 

Spin-charge separation manifests itself within effective 
field theory (bosonization) after linearizing the excitation 
spectrum at the Fermi points. However, the linearization 
brings in additional symmetries, for example Lorentz in- 
variancc, and protects the system from the subtle effect 
of spin-charge interaction. The interaction between spin 
and charge degrees of freedom is seen in exact solutions 
for integrable systems, for instance in the Fermi Hub- 
bard model for spin-1/2 fermions [a, |7|, |8|, |9( . In this case 
the spin-charge coupling can also be treated by bosoniza- 
tion accounting for the curvature of the spectrum at the 
Fermi points [9(. In the presence of two gapless modes, 
this leads to charge transfer by spin excitations |9| . 

One of the achievements of bosonization is the descrip- 
tion of the commensurate-incommensurate (C-IC) phase 
transition for spin-gapped fermions, where the gap is 
closed by a critical magnetic field and a universal square 
root field dep endence of the magnetization emerges at 
half filling [lCj, lll| ■ We show how away from half filling the 
spin-charge interaction leads to a linear critical behavior 
of magnetization and to a finite susceptibility at the C-IC 
phase transition for a fixed number of particles. The spin- 
charge interaction enters the problem through the curva- 
ture of the spectrum at the Fermi points, and our effective 
field theory is applicable for a wide range of models, in- 
cluding continuum and extended Fermi Hubbard models 
with spin anisotropic interactions [12| and/or mass (hop- 
ping) anisotropy [l3|. In fact, for the integrable Fermi 
Hubbard model with only on-site interactions, this type 
of critical behavior was already obtained from the Bethe 
Ansatz [7| . We give a transparent interpretation of this 



picture and show how the spin-charge interaction changes 
the behavior of correlation functions. 

We first consider a dilute strong coupling limit for spin 
gapped fermions and obtain the magnetization across the 
C-IC transition. In this limit spin-| and spin-| fermions 
form strongly bound pairs and the density of the ther- 
modynamic potential is: 
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and d x 9 pi d x 0\ represent density 



The fields d x <fi p , d x 
and current fluctuations for the pairs and uncompen- 
sated (e.g., spin-t) fermions [l|, h is the magnetic field, 
the term W mimics a gap for spin excitations which 
is closed by a critical field h cri [i is the chemical po- 
tential, and the multiple M = (d x (/>i + 2d x <j) p ) / ^/tt de- 
scribes fluctuations of the total number of fermions. At 
a fixed /i the fields d x (f> p , d x 9 p and d x (f>^, d x 9^ are decou- 
pled and one obtains the usual sine-Gordon like square 
root dependence of magnetization on the field [10|, \lM- 
m ~ \/h — h cr , for h — ► h cr + 0. At a constant number 
of particles we have a constraint (J\f) = 0, which pro- 
vides coupling between the fields of spin-f fermions and 
pairs and, as we will see, modifies this square root depen- 
dence to a linear one [14[. At the critical field h cr = 2 A, 
where 2A is equal to the binding energy of the pairs, the 
low-momentum dispersion relation for spin-| fermions is 

£ T (fc) = 
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v 2 k 2 /2A, with v-\ being 



their velocity. The bound pairs disperse linearly with ve- 
locity v p ^ 0. Minimizing il ([T]) with respect to d x <f> p and 
/U and expressing fi through magnetization we obtain: 



2irv\m = \/A(h — h cr — v p Trm). 



(2) 



The appearance of the term linear in m under the square 
root is due to resolving the constraint. Equation @ gives 
m,= (h—h cr ) /■KVp for ft,— >/i cr +0, and the susceptibility is 



X = dm/dh\ hcr = l/irvp. 



(3) 



We now turn to the spin and charge basis and de- 
rive from microscopic principles an asymptotically exact 



theory near the critical point. Taking into account the 
curvature of the spectrum at the Fermi points [9|, Il5| . 
the low-energy Hamiltonian density in the weak coupling 
limit can be written as: 



n = E ^-[(d x <t> a ) 2 /K a + K a (d x e a ) 2 ] 
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+ ^j=k (d x <f> c [(d x <j) s f + (d x 9 s ) 2 ] + 2d x cj)sd x 9sd x 8 c ) , 



where the subscripts c and s stand for the charge and spin 
sectors. The fields d x (j> c and d x 9 c describe fluctuations 
of the charge (mass) density and current, while d x <j) 8 and 
d x 8 s are fluctuations of the spin density and spin current, 
with CjS = Oj ± <j>i)/V2 and C>S = (0 t ± 0{)/\/2. The 
Hamiltonian Q is applicable for a wide range of mod- 
els for spin-1/2 fermions, including continuum and ex- 
tended Hubbard models. The coupling constant g s , Lut- 
tingcr parameters K c ^ s , and deviations of the charge/spin 
velocities v c , s from the Fermi velocity vp depend on 
the Fourier transforms of the interaction potential at 
wavevectors k = and k = 2kp [l|. For spin-gapped 
fermions which arc SU(2) symmetric at h < h cr , one has 
g s < 0, K s = 1 + g s /2, and the charge sector is gapless. 

Compared to the standard bosonized Hamiltonian 
which is quadratic in currents and spin-charge separated, 
Eq.((4]) has an extra (cubic) term [9|, [l5|, [l|| accounting for 
the curvature of the free spectrum at the Fermi points. It 
couples the spin and charge sectors and is proportional 
to k e d 2 E(k)/2dk 2 \k F - The cubic terms within the 
charge sector are omitted in Eq. (0J as they are irrelevant 
modifications of the linearly dispersing charge mode. Di- 
mcnsionally the curvature term is irrelevant as it has the 
scaling dimension equal to 3 > 1 + 1. However, as we 
will argue, this term can not be dropped near the critical 
point where the magnetic field closes the gap, because 
it couples the spin and charge sectors. We will show 
how this coupling effects the magnetic susceptibility at 
the critical point when magnetization sets in. 

For finding the susceptibility at a given number of par- 
ticles we have to impose a constraint: (d x 4> c ) = 0, which 
allows us to integrate out the charge modes. We calcu- 
late the contribution of the curvature term to the ground 
state energy at the onset of magnetization, confining our- 
selves to the terms proportional to m 2 . For extracting 
these terms we write d x cf> s =: d x (j> s : +v27rm, with the 
symbol :: standing for the normal ordering. Then, after 
integrating out charge degrees of freedom, the Euclidean 
action is 5 e // = S® + S K , where S® is the action of the 
sine- Gordon model at magnetization m — > 0, which docs 
not give rise to an m 2 contribution in the ground state 
energy[ 
to to 2 , the term S K originating from the spin-charge in 



teraction is given by: 
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Here r = ivFt is the Euclidean time, x = {x, r} = 
{x ,a;i}, and y = {y,r'} = {y a ,yi\. For simplicity we 
put K c = 1 and v c = v s = vf, which docs not affect our 
main results, and the propagator for the charge sector 
is G c (x,y) = -l/47rln((x - y) 2 /a 2 + (r - r'f/a 2 + 1), 
where a is a short distance cut-off. 

The ground state energy, from which we have to ex- 
tract the m 2 contribution, is given by: 
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where Jdr — ► oo and we have written e~ s " = 1 — S K + 
0(to 4 ). As we see, Eq. © involves the calculation of the 
expectation value of S K §5§ in the vacuum of the sine- 
Gordon theory at m — > 0. 

For h = h cr + the vacuum of effective theory contains 
infinitcsimally small density of solitons. After normal 
ordering with respect to the vacuum at h = h cr + 0, the 
to 2 contribution in Eq. ((5]) can be extracted using the 
vacuum at h = h cr — due to the relation: 

(:d Xi <j> s ::d y:i <f> s :) hcr+0 = (d Xi <t> a d yj <f> s ) hcr _ + O{m). (7) 

Eq. ([7|) can be established from mapping the sine-Gordon 
model onto the massive Thirring model. Then at K s = 
1/2, where the spin sector is equivalent to free mas- 
sive relativistic fermions, one easily gets Eq. J7J). For 
K s y^ 1/2 one finds that Eq. ||7J) holds in any order of 
perturbation theory in the Thirring coupling constant. 

On the other hand, for h < h c the magnetic field does 
not change the states of the system and only shifts the 
antisoliton and soliton energies by ~ ±/i so that the 
energy of a soliton-antisoliton pair remains the same. 
Since d Xi 4> s (yi) has nonzero matrix elements only be- 
tween the states which can differ from each other by a 
certain number of soliton-antisoliton pairs 17j, the cor- 
relation function (d Xi (f> s d V:j cf> s ') for h = h cr — is the 
same as at h = 0. At h = 0, due to the Euclidean in- 
variance we have: \d Xi (f) s dy.<p s ) = d Xi d y .G s (r), where 
r = \J(x — y) 2 + (r — r') 2 is the radial variable. The 
expectation value (S K ) is given by Eq. (O, with the 
products d Xi 4> s {yL)d yi (j) s {y) replaced by the correspond- 
ing correlation functions. Then, after integrating over 
the angular variable, the non-local terms of the first line 
of Eq. ((S]) and those of the second line cancel each other. 
Accordingly, the averaged integrand in Eq. ([5]) reduces 
to: 



[10|,lll|. Retaining only contributions proportional i^ m a Li^'i^^W ^ "^ V -cd 6 (x)d, 
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This means that the m 2 contribution to Eq ((6|) can be 
obtained by using a simplified effective action: 



dxdT[(d x <f> s ) 2 + (d T (f> s ) 2 ].(8) 



S eff : 



One thus sees that the effect boils down to the renor- 
malization of the Luttinger parameter of the spin sector 
(increase of K s ) with m 2 : 



K s -» K s (1 + 2m 2 K 2 TT 2 /v F ) 



(9) 



Note that for spin-gapped fcrmions which are SU(2) sym- 
metric at h < h c Eq. ([9]) encodes breaking of the SU(2) 
symmetry. 

From the rescaling of the Luttinger parameter deter- 
mined by Eq. (|9|) we obtain the following m 2 contribution 
to the ground state energy: 
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For the inverse susceptibility Eq. (flC)]) then yields: 

vlx' 1 = 4K s K 2 n 2 d£ /dK s , (11) 

where So stands for the ground state energy density of the 
sine- Gordon model. For K s — > 1 we can follow the RG 
procedure [18|,[l9| m order to extract the leading universal 
contribution to the ground state energy density of the 
sine-Gordon model. In the one-loop approach we have 
Eq = —XA 2 /vf, where A is the soliton mass (gap), and 
A is a positive factor which we will fix later for the SU(2) 
sine- Gordon case. One- loop RG estimate of soliton mass 
is [H, [H, H : 



A~E F < 
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(12) 



Finally, from Eq. (fTTj) in the vicinity of the SU(2) sepa- 
ratrix we obtain: 
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up to subleading contributions. 

Equation (fl3|) is valid for a wide class of generic mod- 
els, including those with the spin anisotropy. Strictly 
speaking, the Hamiltonian ((4|) requires small g s and K s 
close to unity. Nevertheless, one can think of extending 
our results to K s away from unity, in particular to the 
Luther-Emery point K s — <> 1/2. Then, it is straightfor- 
ward to evaluate x _1 by mapping the spin sector onto 
free massive fermions, which gives: x" 1 oc K 2 d£o/dK s oc 
k 2 ((d t <t> a ) 2 /v 2 F - (d x s ) 2 ) cx ( K Aln A/E F ) 2 . 

The most important result is that the susceptibility 
at the commensurate- incommensurate phase transition 



stays finite if the curvature is finite: k^O. At the onset 
of magnetization the susceptibility is finite also for free 
fermions, where A = and x ~ v f ■ However, as we 
see from Eq. (fl~3|) , in the limit of A — > the susceptibility 
diverges. This was previously observed for the Hubbard 
model [7| and attributed to a singular character of the 
zero interaction point. 

In the case of integrable Fermi Hubbard model with 
only on-site attractive interaction U < 0, one has 1 — 
K s ~ \U\/2ttvf and the result of Eq. (fl~3|) is similar to 
the Bethc Ansatz calculation in the weak coupling limit 
0: X _1 (|^ -*•()) = 8k 2 tt 3 A 2 /v f U 2 . This implies that 
the factor A is equal to 3/27T on the SU(2) line. 

For strong coupling the Bethe Ansatz inverse suscep- 
tibility is given by x~ 1 (\U\ -> oo) = 2tt 2 ^(1 - v) 2 /\U\ 
|7(, which at a low filling factor v tends to our strong 
coupling result ([3]), with v p = 2ttv/\U\. 

We now analyze the behavior of pair and single fermion 
correlation functions at the onset of magnetization (h > 
h cr and m — > 0) [21(. For the Hubbard model, using ex- 
plicitly the dressed charge matrix [22| , in the presence of 
two gapless modes one obtains an effective Hamiltonian 
density @, [23|,[24|]: 



We// = E T [(W 8 /^ + K ^°P) 2 ] 



(14) 



/3=± 



The fields <j>± and 9± are related to the spin and charge 
fields through the spin-charge mixing parameter £: 



<t>+ = <t>c~£,4>s, 9-i 



Os + tfc, (15) 



and v±, K± are the Bethe Ansatz velocities and Lut- 
tinger parameters for the ± sectors. For m — > we have 
V- oc m — > 0, K_ — > 1/2 at any U and v [25j] . In the case 
of half filling (y = 1) one has K + = 1, £ = for all \U\, 
and there is an exact spin-charge separation so that the 
fields <f>±, 9± coincide with <f> c s , 9 CS . For v < 1 one has 



K , = 1 



1^1 
2ttvf 



; ^ 



^cos(^)exp(-^] 



at \U\ — > 0, with the Fermi velocity u^ = 2sin7rz//2 , 
and K+ = 2, f = 1 - i/, v+ = 2nv/\U\ for \U\ -* oo. 
So, the parameter of spin-charge mixing, ^, ranges from 
to (l — i/) and monotonically increases with |[/|. The 
effective Hamiltonian (fT4)) is obtained through the Bethe 
Ansatz calculation and for the inverse susceptibility it 
naturally gives the exact resultQ: X -1 = 27rv + £ 2 ' /K + . 

Asymptotic behavior of correlation functions for the 
Hubbard model with a repulsive on-site interaction, in 
the presence of two gapless modes (and in the presence 
of magnetic field) was obtained by Frahm and Korepin 
[6(. Critical exponents for the general case have been 
obtained from a numerical solution of the coupled Bethe 
Ansatz integral equations for the dressed charge matrix. 
The effective Hamiltonian (TI"4"|) was constructed by Penc 



and Solyom in such a way that it reproduces the Bcthe 
Ansatz behavior of correlation functions [8j . This proce- 
dure of obtaining an effective Hamiltonian was retrans- 
lated to the case of attractive Hubbard model by using 
the particle hole transformation 24 1. 

The limit of rn — ► allows us to derive analytical 
expressions for the critical exponents of the correlation 
functions and make a number of physical conclusions. For 
the pair correlation function from Eq. (fT4|) we obtain: 



(V|(x)^J(x)^(o)Vt(o)) « 



whereas for h < h rr it is 



cos 2-Kmx 

x l/2+l/K + 



oo, (17) 



-i/K + _ There is a universal 



jump of 0.5 in the critical exponent, the result that is ex- 
pected from the theory based on spin-charge separation. 
However, for the single fermion Green function we find: 
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(18) 



where feFfU) is the Fermi momentum of spin- up (down) 
fcrmions given by the free value. The critical exponent of 
the majority (spin-up) component is i/f = 1/2 + A' + /4 + 
(1 + £) 2 /8 + (1 - €) 2 /^K +: and for the spin-down com- 
ponent we obtain v^ = i/f + (1/K + — l/2)£ > v-\. The 
presence of an additional spin-charge mixing term ~ £ 
in the critical exponent of the single fermion Green func- 
tion suggests that u-\ < v± even in the limit of m — > 0, 
which is a clear signature of spin-charge coupling. Per- 
sistence of spin- charge coupling down to m — > limit 
was recently observed numerically [26(. The difference 
v^ — iq increases with \U\ for weak coupling, reaches its 
maximum in the regime of intermediate coupling, and 
then decreases with increasing \U\ in the strong coupling 
regime. Thus, the effect of spin-charge mixing is the most 
pronounced at an intermediate coupling strength. 

In conclusion, we showed that the curvature couples 
spin and charge modes for m — > and changes critical 
properties of ID spin gapped fermions at the onset of 
magnetization. Two-component Fermi gas in a ID op- 
tical lattice is well suited for revealing spin-charge sepa- 
ration or observing spin-charge coupling, especially in a 
box potential where v is coordinate independent. Peri- 
odic modulations of the box size can only excite in-phasc 
oscillations of the two components (charge oscillations), 
and they will not excite out-of-phase oscillations (spin 
mode) at half filling where exact spin-charge separation 
holds. In contrast, for a significantly smaller filling factor 
the excitation of these modes will be provided by spin- 
charge coupling. 
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